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We studied the high driving force regime of the current-voltage transport response in the mixed 
state of amorphous molybdenum-germanium superconducting films to the point where the flux flow 
becomes unstable. The observed nonlinear response conforms with the classic Larkin-Ovchinikov 
picture with a quasiparticle-energy-relaxation rate dominated by the quasiparticle recombination 
process. The measured energy relaxation rate was found to have a magnitude and temperature 
dependence in agreement with theory. 
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I. INTRODUCTION AND THEORY 

In a type II superconductor, magnetizing fields be- 
tween the lower critical field H c \ and upper critical field 
H C 2 introduce flux vortices containing a quantum of flux 
$ = h/2e. A transport electric current density j perpen- 
dicular to the vortices and to exerts a "Lorentz" driving 
force Ft, = j&o- The consequent vortex motion gener- 
ates an electric field E = vB and is opposed by a viscous 
drag r\v (where r\ is the coefficient of viscosity and v the 
vortex velocity), so that in steady state j$> = r]v oc E 
and the response is Ohmic as long as the flow is not 
hindered by pinning. This regime of flux motion corre- 
sponds to free flux flow (FFF). While superficially the 
physics appears simple, and bears resemblance to a hy- 
drodynamic system, this resemblance and apparent sim- 
plicity are deceiving. First of all the so-called "Lorentz" 
force in a superconductor actually has the opposite direc- 
tion to the usual electromagnetic Lorentz force and while 
it has the right magnitude Ft, = j& , the derivation of 
this expression is not completely trivial^—. Second the 
hydrodynamic analogy leads to the expectation that a 
narrower vortex would have a lower drag coefficient. Un- 
der normal circumstances the reality is just the opposite: 
dissipation in and around the core of a vortex arises from 
generated electric fields acting upon normal quasiparti- 
cles leading to Ohmic dissipation^ and from irreversible 
entropy transfer occurring between the leading and trail- 
ing edges of moving vortices^. These dissipativc processes 
increase with the order-parameter gradient so that 7/ is 
roughly inversely proportional to the vortex core area. 
Thus under normal circumstances r\ drops as the the vor- 
tex expands and hence the flux-flow resistance goes up 
with temperature^^. 

The Ohmic FFF regime discussed above ceases as j and 
E are increased to high values that alter the supercon- 
ducting state. If the electron-electron scattering time r ee 
is short compared to the the electron-phonon scattering 
time r ep , the distribution function of the quasiparticles 
becomes thermalized and a finite power dissipation den- 
sity jE simply raises the electron temperature to a value 
T' above the phonon and substrate temperatures, T p and 
Tq respectively This causes an expansion of the vortex 



core and a drop in viscosity leading to a non-monotonic 
("N" shaped) j(E) curve and a consequent instability. 
Such a hot-electron instability model was developed and 
experimentally verified by us in our earlier work^i. 

A more intricate scenario, proposed by Larkin and 
Ovchinniko v 12 ' 13 (LO), arises when r ee >r ep . In this 
case the quasiparticle distribution function acquires a 
non-thermal shape and the vortex dynamics are altered 
in a complicated and less obvious way. At high E the 
quasiparticle population in the core reduces while it in- 
creases outside the core, thus causing the vortex to shrink 
while simultaneously reducing the contrast in quasipar- 
ticle density between the two regions. Thus the gradient 
in the order parameter A is not boosted to the extent an- 
ticipated by the reduced vortex size. In the meanwhile, 
the diminished quasiparticle population in the core tends 
to lessen the dissipation and reduce the drag. Thus over- 
all rj declines despite the shrinking in size with increasing 
E. Like the hot-electron case discussed earlier, this again 
leads to a non-monotonic j(E) curve and a vortex insta- 
bility. The LO instability has been observed in previous 
experiments^— and the combination of heating effects 
and the LO mechanism have been considered by various 
authors 2 ^— . 

The relative magnitudes of r ee and r ep govern 
which mechanism dominates the instability. Standard 
estimates 3 ^ for the scattering times r ee =rhtF /k B T 2 and 
T ep =r 3 hTp/kBT 3 (where Try is the Debye temperature 
and r < 1 is the phonon reflection coefficient at the film- 
substrate interface arising from acoustic mismatch) give 
a cross-over temperature of Tx = r 2 kBTp/ep. In the 
present work, we investigated amorphous molybdenum- 
germanium (MoGe) films, which have (estimating from 
known parameter a 31 ' 32 ) Td~260 K and €f ~ 10 cV, and 
thus have T x < 0.6r 2 K. Since r < 1, T x will be well 
below the temperature range of our experiment (3-6 K). 
Thus we expect the non-linear j(E) response to be dic- 
tated by the LO mechanism, which is indeed born out by 
our data. 

In the LO theory, the non-linear flux-flow conductivity 
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at high E is given by (Eqs. 38 and 53 of Rcf. [l 



a(E) » ^ 
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in terms of the free- flux- flow value <r/ of the linear regime 
(_E — > limit in the absence of pinning) and E* the crit- 
ical electric field at which j attains its maximum value 
before the j(E) curve enters a region of negative differ- 
ential conductivity, t = T/T C is the reduced temperature 
(T c is the superconducting transition temperature), c is 
an unknown constant of order unity, and E* is given by 
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where D is the diffusion constant, Q{x) is the Riemann 
zeta function, and r e is the energy relaxation time. 

In our previous work^ we found that the expressions 
for cr/(0) in reference UH (their Eqs. 22 and 30) did not 
fit the data well over a significant range. Instead the 
following expression^ 
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based on the mean-field result of time-dependent 
Ginzburg Landau (TDGL) theory^—, more accurately 
represented the behavior over an extended range; here 
b = B/hqH C 2 is the reduced magnetic field (H c2 is the 
upper critical magnetizing field) and v ~ 0.3 is a dimen- 
sionless constant. The right hand side of Eq. [3] represents 
a two-fluid-model sum of the normal conductivity a n and 
the flux-flow conductivity contribution. The first term is 
negligible compared to the second term for the range of 
conditions where we study the LO nonlinear effect (i.e., 
our mixed-state conductivity is far higher than the nor- 
mal conductivity) and also the a n term remains constant 
and is not affected by E and B fields. Thus combining 
Eq. [1] with the second term of Eq. [3] gives the following 
nonlinear j(E) relationship 



(4) 



Note that this non-linear function is only valid until the 
vortex stops shrinking, which occurs at a fielcU^ E s ~ 
E*/(l — t) 1 / 4 . At very high E the system eventually 
enters the normal state and then a = a n . 



II. EXPERIMENTAL TECHNIQUES 

The samples A, B, and C used in this experiment are 
exactly the same as the samples A, B, and C used in 
our prior work on free flux flow^. The samples consist 
of MoGe films of thickness 50 nm sputtered onto sili- 
con substrates with 200 nm thick oxide layers using an 




alloy target of atomic composition Moo.79Geo.21- The 
deposition system had a base pressure of 2 x 10~ 7 Torr 
and the argon-gas working pressure was maintained at 3 
mTorr during the sputtering. The growth rate was 0.15 
nm/s. The samples were patterned into bridges of length 
/ = 102 /im and width w — 6 /im using photolithog- 
raphy and argon ion milling. Some parameters of the 
samples are as follows: Sample A: T c =5.56 K, i?„=555 
Q, (i H' c2 =-3.13 T/K and D = 0.35 cm 2 /s. Sample B: 
T c =5.41 K, i?„=555 ft, n H' c2 =-3.13 T/K and D = 0.35 
cm 2 /s. Sample C: T c =5.01 K, i?„=630 Q, fi o H' c2 =-3.0 
T/K and D = 0.37 cm 2 /s. Here, R n is the normal- 
state resistance, H' c2 = dH c2 /dT\r c is the upper-critical- 
field slope, and the diffusion coefficient D was calculated 
fromia D = ~8k B /2ne^ H' c2 . 

The cryostat was a Cryomech PT405 pulscd-tube 
closed-cycle refrigerator that went down to about 3.2 
K. It was fitted inside a 1.3 Tcsla GMW 3475-50 water- 
cooled copper electromagnet. Calibrated cernox and hall 
sensors monitored T and B respectively. The electri- 
cal transport measurements were made with an in-housc 
built pulsed current source, preamplifier circuitry, and a 
LeCroy model 9314A digital storage oscilloscope. The 
pulse durations are 20 /is or less, with duty cycles in the 
few ppm range to reduce macroscopic heating of the film. 
Our previous review paper a 41 ' 42 give further details about 
the measurement technique and the thermal analysis. 



III. DATA AND ANALYSIS 

Fig. [T] shows some examples of nonlinear j(E) curves. 
j (and hence the viscous drag force) has a local maxim- 
imum value of j* at the instability field E* . In a cur- 
rent biased circuit, where the source resistance is larger 
than the sample resistance as is the case here, E jumps 
(indicated by arrows) upon increasing j to the vicinity 
of the intrinsic j*. In a voltage biased measurement, 
where the source resistance is lower than the sample's, 
the there will not be a jump in E and instead j will 
be seen to decrease. The macroscopically averaged be- 
havior will have a negative dj /dE and the flux matter 
fragments into compressional 4 ^, shear— or other types 
of elastic domains such that any given domain is mov- 
ing in a response region with dj /dE > locally. The 
macroscopic j(E) curve will then not follow the primi- 
tive curve (e.g., Eq. [J) but will show steps in the region 
where dj /dE < 0. In the present experiment we are only 
concerned with the region of the transport response up 
to E* where dj/dE > macroscopically. 

The solid lines in Fig. [T] arc fits to Eq.|4]and are seen to 
follow the trends of the data. The parameters v and E* 
(location of peak) were adjusted to improve the fits, but 
have fit values of the expected magnitudes: v ~ 0.3 and 
E* from the peaks is slightly higher than the position of 
the actual jumps, which is to be expected and has been 
observed by others (e.g., reference \2u) . For subsequent 
analysis, we take E* to be the actual measured value of 
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E at the threshold of the jump. 

Fig. |H plots E* 2 , obtained from the j(E) curves as 
discussed above, against B 2 . In agreement with Eq. [21 
the two quantities are directly proportional to each other 
(i.e., the critical vortex velocity v*=E* /B is independent 
of B). From the measured slope and Eq.fJJ we obtain the 
energy relaxation time r s . Fig.[3]plots the corresponding 
relaxation rate r" 1 against the reduced temperature for 
each of the three samples. 
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FIG. 1: j versus E curves for sample B at two temperatures 
(a) 3.40 and (b) 4.63 K, at indicated B values. The symbols 
correspond to the measured data and the solid lines corre- 
spond to Eq. [4] The dashed arrows indicate the observed 
jump in E. For the theoretical curves, v and E* were ad- 
justed to the following values: For T =3.40 K, v ~ 0.3 and 
E*= 34, 68, 138, and 292 V/m in the order of ascending B. 
For T=4.63 K, v » 0.2; and E*= 20, 40, 83, and 171 V/m. 

As discussed in the introduction, the LO effect occurs 
when T ee is long compared with r ep , resulting in a non- 
thermal shape of the quasiparticle distribution function. 
The extent of the distribution function distortion is con- 
trolled by the rate of energy relaxation from quasiparti- 
cles to phonons, which occurs through two processes: one 
process is the inelastic scattering between a quasiparticle 
and a phonon and the other is the recombination of two 
quasiparticles to form a Cooper pair with the emission of 
a phonon. As discussed by Kaplan et al4^, the energy re- 
laxation is mainly dominated by the latter recombination 
process which has a rate that can be written as: 
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FIG. 2: E* 2 versus B 2 for samples A, B, and C at indicated 
temperatures showing proportionality between the two quan- 
tities as per Eq. [2] The lines are least-squares linear fits to 
the data (symbols). 



where T is a temperature independent characteristic 
time constant (in the terminology of Kaplan et al4^, 
7" ~ To/55) and A is the temperature dependent su- 
perconducting energy gap. Taking the BCS (Bardeen- 
Cooper-Schrieffer) temperature dependence for A we are 
able to fit the data in Fig. [3] with Eq. [5] (solid lines) with 
T ~ 3.5 x 10~ n s. As can be seen, the measured func- 
tional form of r £ _1 (t) is in agreement with Eq. [5j While 
there is insufficient information in the literature to theo- 
retically compute the magnitude of T for comparison, Ta- 
ble I of Reference lists values of To (~ 55 x T) for vari- 
ous other materials. Although there is an enormous range 
in T — from 8xl0~ 13 to4xl0~ 5 sfor various materials — 
it is interesting to observe that tantalum with about the 
same T c (4.5 K) and T D (240 K) as MoGe (for which 
T c «5.3 K and T D «260 K) has a value of T (3.3 x 1CT 11 
s) that is comparable to the one we obtained for MoGe 
(3.5 x 10 -11 s). (T c and Tn are parameters that are 
indicative of the electron-phonon coupling and phonon 
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density of states respectively.) If quasiparticlc-phonon 
scattering is the dominant relaxation process, rather than 
quasiparticlc recombination, then the rate is given by a 
function of the formic r^ 1 oc i 7 / 2 instead of Eq. [5] This 
power-law function (taken with an adjustable constant 
of proportionality a) is plotted as dashed red lines on 
Fig. |3] and is clearly at odds with the data. Thus our 
study of the vortex instability is able to distinguish the 
two routes of energy decay and provides a confirmation 
of the recombination rate expression of Eq. [5j 
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FIG. 3: Temperature dependence of the energy-relaxation 
rate. The symbols represent the experimental values deduced 
from Eq. [2] and the measured E* . The solid black lines are 
fits to Eq. [S] (quasiparticle recombination process) with the 
values T = 3.4 x KT 11 s, 3.7 x KT 11 s, and 3.3 x KT 11 
for samples A-C respectively. The dashed red lines are fits 
to r^ 1 = at 7 ^ 2 (expected for inelastic quasiparticle-phonon 
scattering^) with a = 8 x 10 9 s"\ 7 x 10 9 s" 1 , and 8 x 10 9 
s _1 for samples A-C respectively. 



IV. CONCLUSIONS 



In conclusion, we have studied the high driving force 
regime of vortex dynamics in one of the simplest and 
nearly model superconductors (unpinned, isotropic, low- 
temperature, weak-coupling BCS, etc.). In recent work2£ 
we found that for the FFF regime these MoGe films pro- 
vided a detailed confirmation of the TDGL mean-field 
prediction for 07 (0), while the LO expressions for the 
same regime showed very limited applicability. On the 
other hand in the present work we find that the LO ex- 
pression (Eq. [1]) for the nonlinear modulation factor for 
tr/(0) and the LO result (Eq. [5J for the relationship be- 
tween the instability electric field E* and the energy re- 
laxation time r e , are well obeyed in this system. Further- 
more we were able to distinguish between the two princi- 
pal quasiparticle-phonon energy relaxation processes and 
confirm the predicted temperature dependence (Eq. [5]) 
for the recombination process. Both regimes of instabil- 
ity, the hot-electron as well as the distribution-function 
type, thus provide a valuable tool for investigating key 
scattering processes through the response of the mixed 
state. 
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